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ABSTRACT 

An  interactive  program  has  been  installed  at  the  Fort  Collins  Computer 
Center  (FCCC)  to  assist  users  to  create  data  files  required  for  the 
Douglas-fir  Tussock  Moth  Stand.  Outbreak  Mod.el  programs.  Instructions  on  how 
to  use  this  interactive  program  are  provided  via  examples  in  this  report. 


INTRODUCTION 

Through  a cooperative  effort  between  Pacific  Northwest  Forest  and  Range 
Experiment  Station  and  FPM/Methods  Application  Group,  a report  entitled  "Data 
Preparation  and  Computer  Runstream  Procedures  for  the  Douglas-fir  Tussock  Moth 
Stand  Outbreak  Model"  was  distributed  (Colbert  and  Wong  1979).  It  was 
prepared  under  the  premise  that  a comprehensive  reference  document  should  be 
assembled  to  serve  as  a single  source  of  information  on  data  input  and  output 
structure  and  computer  access  procedures  for  this  model.  The  stand  outbreak 
model  is  a component  of  the  Douglas-fir  Tussock  Moth  (DFTM)  Pest  Management 
System  (Campbell  and  McFadden)  developed  under  the ‘USDA  Expanded  DFTM  Research 
and  Development  Program  (Wright  1977).  This  model  has  a generalized  system 
structure.  When  new  information  becomes  available  it  can  be  readily  updated. 

The  basic  documentation  on  the  DFTM  Stand  Outbreak  Model  described  data 
input  and  output  structure  for  a variety  of  available  options.  However,  the 
need  exists  to  further  simplify  the  data  preparation  procedures  so  that  users 
with  a minimal  amount  of  training  can  readily  create  data  files  required  to 
make  simulations.  To  satisfy  this  need,  an  i nteracti ve  program  has  been 
developed  and  placed  on  line  at  FCCC.  Designed  strictly  for  demand 
processing,  this  program  can  be  accessed  by  an  authorized  FCCC  user  from  a 
remote  terminal . 

In  this  report,  it  is  assumed  that  the  reader  has  already  reviewed  the 
aforementioned  documentation  for  the  DFTM  Stand  Outbreak  Model  (Colbert  and 
Wong  1979).  The  reader  should  therefore  be  knowledgeable  about  the  input 
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files  required  and  the  information  needed  to  make  up  these  files.  Primary 
intent  of  this  report  is  to  demonstrate  how  this  new  program  can  be  utilized. 
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SYSTEMS  OVERVIEW 

The  interactive  program  is  conversational  in  nature  and  is  capable  of 
performing  the  following  tasks: 

1.  Create  the  initial  condition  file  (*IC.)  and  the  parameters  file 
(*PARAMETERS. ) simultaneously  (option  1)  as  input  files  for  a 
simul ation . 

2.  Create  the  specifications  file  {*SPECTAR4.)  to  retrieve  detail 
information  from  a previous  simulation  (option  2). 

3.  Compute  instar  specific  daily  larval,  pupal  or  overwintering 
mortality  resulting  from  direct  control  as  input  to  the  generation  of 
the  parameters  file  (option  3). 

This  program  is  written  in  ASCII  Fortran  for  the  UNIVAC  1100  computer  at 
FCCC.  It  resides  in  a program  file  named  DFTM*UTILITY1 . 


PROGRAM  FEATURES 


Several  special  features  have  been  built  into  the  interactive  program 
v/hich  will  transform  field  data  into  a form  readily  acceptable  by  the  model 
programs.  Some  of  these  special  features  are: 

1.  Capability  to  Enter  Model  at  Any  Stage  of  the  Outbreak. 

The  DFTM  Stand  Outbreak  Model  was  structured  in  a form  such  that  the 
insect  density  (larval  count)  and  the  actual  foliage  conditions  must  be 
specified  at  the  initiation  (1st  in  star)  of  a phase.  This  requirement  is 
not  practical  from  the  standpoint  that  insect  densities  at  the  1st  instar 
are  not  alv/ays  available  except  under  an  established  monitoring  system 
when  field  activities  are  planned  in  advance.  Recognizing  this  situation, 
the  capability  to  "back-up"  insect  count  available  at  some  intermediate 
point,  e.g.  at  the  third  day  of  the  second  instar,  to  the  beginning  of  the 
first  instar,  has  been  incorporated  into  this  program.  To  accomplish  this 
task,  the  following  functional  form  for  insect  survivorship  is  utilized: 

fr  fr 

i = l 


j = l 


(1  - 
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where  N = insect  density  at  the  initiation  of  a specific  phase. 
Nq  = instar  density  at  the  time  of  data  collection, 
i = the  type  of  mortality  associated  with  instar  i. 

P-jj  = the  mortality  (value)  associated  with  i and  j. 
d-jj  = the  time  in  days  to  be  included  for  each  instar. 

2.  Ability  to  Adjust  Mortality  Factors  and  Foliage  Biomass  Values. 


A set  of  natural  mortalities  are  stored  within  the  program  in  the 
*PARAMETERS.  file.  Users  can  utilize  these  standard  mortalities  or  make 
alterations  when  executing  this  program  if  other  values  are  to  be  used. 
Foliage  condi tions--actual  new  foliage  biomass  and  actual  old  foliage 
biomass  in  addition  to  nominal  percent  new  foliage  and  nominal  total 
foliage  biomass--at  the  first  instar  of  a specific  phase  must  be  available 
for  a simulation  to  start  from  other  than  phase  1.  For  the  phase  1 start, 
the  actual  foliage  complements  will  be  automatically  computed  from  the 
nominal  foliage  values  supplied  by  the  user. 

3.  Ability  to  Create  an  Input  File  to  Retrieve  Detailed  Information. 

Option  2 provides  the  mechanism  to  set  up  the  SPECTAB4.  file  to  retrieve 
model  details  for  up  to  six  state  or  output  variables;  namely,  new  foliage 
biomass,  old  foliage  biomass,  number  of  insects,  mean  larval  biomass, 
number  of  days  on  old  foliage,  and  number  of  days  without  food.  By 
entering  the  appropriate  pair  of  indices  sequentially,  as  shown  in  Fig.  8 
of  the  basic  considerations,  the  desired  output  can  be  obtained  by 
executing  DFTM*RUNSTREAM.TABLE4B. 

4.  Control  Mortal ity  Computation. 

Another  special  feature  is  found  in  option  3,  which  is  the  control 
mortality  algorithm.  With  this  option,  the  user  can  obtain  a value  for 
the  required  control  mortality  factor  to  be  entered  into  *PARAMETERS,  when 
option  1 is  processed.  Using  the  same  functional  form  as  in  the  above 
formula,  the  daily  control  mortality  or  overwintering  and  pupal  control 
mortality,  as  applicable,  can  be  easily  calculated.  For  the  daily  control 
rate,  it  is  assumed  that  the  initial  and  final  insect  densities  are  known 
at  the  beginning  and  end  of  instar.  It  is  then  computed  as  a function  of 
these  values,  and  all  other  mortalities  entered.  Furthermore,  if  stress 
mortality  is  to  be  considered,  in  computing  daily  control  mortality,  the 
user  is  expected  to  supply  the  number  of  days  for  which  stress  will 
operate^.  In  order  to  obtain  an  estimate  of  the  number  of  days  for 
which  stress  will  operate  a preliminary  simulation  will  have  to  be  made 


3 Number  of  days  stress  operates  is  the  number  of  days  on  old  foliage. 
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v/ith  a phase  1 start,  with  insect  densities  ranging  from  1-20  larvae/1000 
square  inch  of  foliage,  and  with  nominal  foliage  conditions  that  best 
approximate  the  situations  at  hand.  After  this  simulation  is  made,  the 
information  on  the  number  of  days  on  old  foliage  can  be  retrieved  as 
discussed  above. 


EXAMPLES 


The  following  synopsis  was  developed  to  show  an  example  of  how  the 
Douglas-fir  Tussock  Moth  Stand  Outbreak  Model  with  its  new  interactive  front- 
end  program  could  be  used  to  make  a pest  management  decision. 

A combination  of  pheromone  trapping  and  larval  sampling  indicates  that  a 
major  up-surge  of  Douglas-fir  tussock  moth  is  about  to  occur  in  portions  of 
eastern  Oregon  during  late  1983.  The  infested  area  encompasses  approximately 
150,000  acres  of  predominately  grand  fir  type.  The  outbreak  is  assumed  to  be 
entering  Phase  II  in  1984.  In  October  1983,  a decision  must  be  made  whether 
or  not  to  begin  plans  for  a direct  control  project. 

Data  from  biological  evaluations  identifies  five  areas  of  infestation. 
Projected  average  1st  instar  population  trends  in  these  areas  at  the  beginning 
of  Phase  II  are  as  follows: 


Infestation  Area 


Projected  1st  Instar  Larvae 
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At  the  time  of  this  initial  evaluation,  no  data  are  available  on  incidence 
of  naturally  occuring  virus  in  the  population  of  egg  parasite  levels  or 
porportion  of  nonviable  eggs.  According  to  existing  data,  the  proportion  of 
nominal  new  foliage  for  grand  fir  in  eastern  Oregon  is  30  percent  and  nominal 
foliage  biomass  on  a mid-crown  branch  is  230  grams/1000  sq.  in.  The  actual 
amount  of  new  and  old  foliage  at  the  beginning  of  Phase  II  is  estimated  to  be 
69  grams  and  161  grams,  respectively. 

One  alternative  under  consideration  is  direct  control  of  the  infestation 
with  an  aerial  application  of  Sevin  4-Oil  at  2 lbs  to  3/4  gal.  carrier  per 
acre.  Data  from  a 1974  pilot  control  project  in  northern  Idaho  (Ciesla  et  al . 
1976)  is  used  to  compute  daily  mortality  rates;  and  the  Stand  Outbreak  Model 
is  run  with  or  without  introduction  on  direct  control  to  project  population 
trends  and  impact.  Data  from  this  pilot  control  project  were  collected  under 
the  following  conditions: 
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1.  Phase  of  outbreak  is  Phase  III. 

2.  Instar  to  v/hich  spray  was  applied  was  primarily  the  3rd 
instar. 

3.  Prespray  population  density  was  50.71  insects  per  1,000  sq. 
in.  of  fol iage. 

4.  Seven-day  post-spray  population  density  = 5.58  insects  per 
1,000  in.  of  foliage. 

5.  Fourteen  day  post-spray  population  density  = .96  insects  per 
1,000  in.  of  fol iage. 

To  use  the  Stand  Outbreak  Model  for  pest  management  decisionmaking, 
compute  a daily  control  mortality  rate  from  the  1974  data;  enter  that  rate 
into  the  parameter  file;  and  run  simulations  with  and  without  a direct  control 
mortality  factor  with  a Phase  II  start;  and  then  compare  the  population  trends 
and  impact. 

The  first  step  is  to  create  *IC.  and  *PARAMETERS.  files  for  the  no  control 
option  as  follows: 
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Runstream  No.  1 


Using  these  input  files,  make  a simulation  without  control  as  follows: 

>0RUN  M01MAG, 11052034050  ,W0NG 
DATE:  080680  TIME:  114337 

S2K  INFO  - , $2000,80=00006  (@I NF0)M0NDAY  11:35 

<(3XQT  DFTM*UTILITY1. 

DFTM  OUTBREAK  MODEL  DATA  PREPARATION  UTILITY 
VERSION  01.001  COMPILED  6/3/80 


PROGRAM  OPTIONS:  1=  IC  AND  PARAMETERS  FILE  GENERATION 
2=  SPECTAB4  FILE  GENERATION, 

3=  CONTROL  MORTALITY  ALGORITHM 
ENTER  OPTION  DESIRED 
>1 

*IC.  FILE  AND  *PARAMETER.  FILE  GENERATION 

ENTER  RUN  QUALIFIER 

>WC1 

RUN  QUALIFIER  IS  MOW  WCl 
ENTER  NUMBER  OF  DOUGLAS  FIR  TREE  CLASSES 
>0 

ENTER  NUMBER  OF  GRAND  FIR  TREE  CLASSES 
>5 

ENTER  STARTING  PHASE  NUMBER 
>2 

DO  YOU  WANT  TO  DEFINE  SUBSETS?  (YES/NO) 

>N0 

DO  YOU  WISH  TO  CHANGE  PARAMETERS  FILE  VALUES?(YES/MO) 
>N0 


INITIAL  CONDITIONS  FILE  GENERATION 
ENTER  NUMBER  OF  DAYS  STRESS  OPERATES  IN 
FIRST  INSTAR,  SECOND  INSTAR,  THIRD  INSTAR,  FOURTH  INSTAR 
>0,0, 0,0 


TREE  CLASS  DEFINITIONS: 

YOU  HAVE  SPECIFIED  000  DOUGLAS  FIR  TREE  CLASSES 
AND  005  GRAND  FIR  TREE  CLASSES.  FOR  EACH  CLASS 
YOU  WILL  NEED  TO  ENTER  1 LINE  CONTAINING  THE  FOLLOWING  VALUES: 

1.  NOMINAL  % NEW  FOLIAGE 

2.  NOMINAL  TOTAL  FOLIAGE  BIOMASS 

3.  # OF  VIABLE  EGGS  PER  BRANCH 

4.  INSTAR  OF  EGG  COUNT 

5.  DAY  IN  INSTAR  OF  EGG  COUNT 
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6.  ACTUAL  NEW  FOLIAGE  BIOMASS 

7.  ACTUAL  OLD  FOLIAGE  BIOMASS 


>30,230,50,1,1,69,160 

>30,230,75,1,1,69,161 


>30,230,110,1,1,69,160 

>30,230,150,1,1,69,161 

>30,230,200,1,1,69,161 

IF  YOU  WANT  A LISTING  OF  THE  ENTIRE  *IC  FILE,  ENTER  -1,-1,0 

IF  YOU  WANT  TO  SEE  ONE  LINE,  ENTER  LINE#, 0,0 


IF  YOU  WANT  TO  CHANGE  A VALUE,  ENTER  LINE#, ITEM#, NEW  VALUE 
TO  CONTINUE  WITH  PROGRAM,  ENTER  0,0,0 


>1,7,161 
LINE  1 ITEM 
>3,7,161 
LINE  3 ITEM 
>-l,-l,0 
2 .00000 
2 .00000 
2 .00000 
2 .00000 
2 .00000 
>0,0,0 


7=  161.000 
7=  161.000 

.00000  30.00000230.00000 
.00000  30.00000230.00000 
.00000  30.00000230.00000 
.00000  30.00000230.00000 
.00000  30.00000230.00000 


69.00000161.00000  50.00000 

69.00000161.00000  75.00000 

69.00000161.00000110.00000 
6 9 . 00000 161. 000001 50 . 00000 

69 . 000001 61 . 00000200 . 00000 


*IC.  AND  *PARAMETERS.  FILE  PREPARATION  FOR  PHASE  2 START  COMPLETED 


****QQQ[)g 

>0FIN  4 


1 

2 

3 

4 

5 


4 Cost  for  this  run  was  $1.34 
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Using  these  input  files,  make  a simulation  without  control  as  follows: 


Runstream  No.  2 


>mm  M01MAG,  1105203405  , WONG 

DATE:  080680  TIME:  115321 

S2K  INFO  - , 52000,80-00006  (01 NF0)M0NDAY  11:35 

>0QUAL  WCl 

READY 

>0ADD,L  DFTM*RUNSTREAM.PRE2 


0DELETE,C  *LUN25. 

FURPUR  27R3A  E35  SL73R1  08/06/80  11:54:34 

*LUN25  IS  NOT  CATALOGUED  OR  ASSIGNED 

FAC  STATUS:  400010000000 


0ASG,UP  *LUN25. 
READY 


*USE  25,*LUN25 
READY 


0DELETE,C  *LUN10. 

FURPUR  27R3A  E35  SL73R1  08/06/80  11:54:44 

*LUN10  IS  NOT  CATALOGUED  OR  ASSIGNED 
FAC  STATUS:  400010000000 


OASG,UP  *LUN10. 
READY 


OUSE  1O,*LUN10 
READY 


OASG,T  *4 
READY 


OASG,T  *5 
READY 


OASG,T  *6 
READY 


PASG,!  *7 
READY 


0ASG,T  *8 
READY 


0ASG,T  *9 
READY 


0ASG,T  *11 
READY 


@ASG,T  *12 
READY 


0DELETE,C  *15. 

FURPUR  27R3A  E35  SL73R1  08/(!!6/80  11:55:34 

*15  IS  NOT  CATALOGUED  OR  ASSIGNED 

FAC  STATUS:  400010000000 


0ASG,UP  *15 
READY 


0FREE  *15. 
READY 


0ASG,A  *15. 
READY 


0DELETE,C  *16. 

FURPUR  27R3A  E35  SL73R1  08/06/80  11:56:06 

*16  IS  NOT  CATALOGUED  OR  ASSIGNED 

FAC  STATUS:  400010000000 


0ASG,UP  *16 
READY 


0FREE  *16 
READY 


*ASG,A  *16 
READY 
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*DELETE,C  *17 

FURPlJR  27R3A  E35  SL73R1  08/06/80  11:56:28 


*17  IS  NOT  CATALOGUED  OR  ASSIGNED 

FAC  STATUS:  400010000000 


0ASG,UP  *17 
READY 


0FREE  *17. 
READY 


PASG,A  *17. 
READY 


0ASG,T  *18 
READY 


0ASG,T  *19 
READY 


0ASG,T  *20 
READY 


0DELETE,C  *21. 

FURPUR  27R3A  E35  SL73R1  08/06/80  11:56:52 

*21  IS  NOT  CATALOGUED  OR  ASSIGNED 
FAC  STATUS:  40010000000 


0ASG,UP  *21 
READY 


0FREE  *21. 
READY 


0ASG,A  *21 
READY 
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0ASG,T  *22 
READY 


(3ASG,T  *23 
READY 


0ASG,T  *24 
READY 


@ASG,A  DFTM*DFTMDF. 

FACILITY  WARNING  000200000000 


0USE  13,DF'm*DFTMDF 
READY 


0ASG,A  DFTM*DFTMGF. 

FACILITY  WARNING  000200000000 


0USE  14,DFTM*DFTMGF 
READY 


0ASG,A  *IC. 
READY 


0ASG,A  *PARAMETERS. 
READY 


0USE  30,*IC 
READY 


0USE  60,*PARAMETERS 
READY 

>0XQT  DFTM*OUTBREAK. 

>0ADD,L  DFTM*RUNSTREAM.P0ST2 


0FREE  *2 

FACILITY  WARNING  100000000000 


@FREE  *3 

FACILITY  WARNING  100000000000 


0FREE  *4 
READY 


(BFREE  *5 
READY 


0FREE  *6 
READY 


0FREE  *7 
READY 


0FREE  *8 
READY 


PFREE  *9 
READY 


gFREE  *10 
READY 


PFREE  *11 
READY 


0FREE  *12 
READY 


PFREE  *13 
READY 


PFREE  *14 
READY 


PFREE  *18 
READY 


PFREE  *19 
READY 


PFREE  *20 
READY 
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@FREE  *21 
READY 


0FREE  *22 
READY 


0FREE  *23 
READY 


@FREE  *24 
READY 


0FREE  *25 
READY 

>0SYM,U  *15.,,FCR104 
>0SYM,U  *16.,,FCR104 
>0SYM.U  *17.,,FCR104 

>0FIn5 


5 Cost  for  this  run  was  $3.37. 
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Then  compute  a direct  control  daily  mortality  rate.  The  runstream  appears  as 
fol  1 ows : 


Runstream  No.  3 


UMSX9RUM  M01KAG,  1105203405  WONG 

DATE:  080680  TIME:  113813 

S2K  INFO  52000,80-00006  ( 0INF0)M0NDAY  11:35 

>0XQT  DFTM*UTILITY1. 

DFTM  OUTBREAK  MODEL  DATA  PREPARATION  UTILITY 
VERSION  01.001  COMPILED  6/3/80 


PROGRAM  OPTIONS:  1=  IC  AND  PARAMETERS  FILE  GENERATION 

2=  SPECTAB4  FILE  GENERATION, 

3=  CONTROL  MORTALITY  ALGORITHM 
ENTER  OPTION  DESIRED 
>3 

CONTROL  MORTALITY  ALGORITHM 

ENTER  OPTION  CODE 

DAY=DAILY  CONTROL  MORTALITY 
PUP=PUPAL  CONTROL  MORTALITY 
OVW=OVERWINTER  CONTROL  MORTALITY 
TER=TERMINATE  PROGRAM 

>DAY 

ENTER  PHASE 
>3 

ENTER  INSTAR 
>3 

ENTER  TREE  TYPE (DF=DOUGLAS  FIR,GF=GRAND  FIR) 

>GF 

NATURAL  DAILY  MORTALITIES: 

PREDATOR/PARASITE=  .003  DISEASE=  .006  BACK  GROUND=  .020  STRESS=  .100 
CHANGE  ANY  NATURAL  MORTALITIES? ( YES/NO) 

>MO 

ENTER  NUMBER  OF  DAYS  STRESS  OPERATES 

>0 

ENTER  BEGINNING  INSECT  COUNT 
>50.71 

ENDING  INSECT  COUNT  WITH  NATURAL  MORTALITY  ONLY=  37.9 

ENTER  ENDING  COUNT  WITH  CONTROL 

>3.06 

DAILY  CONTROL  MORTALITY=  .2239 

ENTER  OPTION 

>TER 

ROUTINE  TERMINATED 

CONTROL  MORTALITY  ALGORITHM  COMPLETED 
****GQODBYE**** 

>OFIN/ 


6 This  number  was  estimated  from  1974  pilot  project  data 

7 Cost  for  this  run  was  $0.37. 


Next  create  the  *IC.  and  *PARAMETERS.  File  for  the  direct  control 
option  as  follows: 


Runstream  No.  4 


>0RUN  M01MAG, 1105203405  ,W0NG 
DUP  ID,  NEW  ID  IS  M01MAH 
DATE:  080680  TIME:  120321 

S2K  INFO  - , 52000,80-00006  (01 NF0)M0NDAY  11:35 
>0XQT  DFTM*UTILITY1. 


DFTM  OUTBREAK  MODEL  DATA  PREPARATION  UTILITY 
VERSION  01.001  COMPILED  6/3/80 


PROGRAM  OPTIONS:  1=  IC  AND  PARAMETERS  FILE  GENERATION 
2=  SPECTAB4  FILE  GENERATION, 

3=  CONTROL  MORTALITY  ALGORITHM 
ENTER  OPTION  DESIRED 
>1 

*IC.  FILE  AND  *PARAMETER.  FILE  GENERATION 

ENTER  RUN  QUALIFIER 

>WC2 

RUN  QUALIFIER  IS  NOT  WC2 
ENTER  NUMBER  OF  DOUGLAS  FIR  TRE  CLASSES 
>0 

ENTER  NUMBER  OF  GRAND  FIR  TREE  CLASSES 
>5 

ENTER  STARTING  PHASE  NUMBER 
>2 

DO  YOUT  WANT  TO  DEFINE  SUBSETS?  (YES/NO) 

>NO 

WHAT? 

>NO 

DO  YOU  WISH  TO  CHANGE  PARAMETERS  FILE  VALUES?(YES/NO) 

>YES 

TO  CHANGE  A PARAMETERS  FILE  VALUE, ENTER  A TRIPLE 
(RECORD  NUMBER, ITEM  NUMBER, NEW  VALUE) 

TO  OBTAIN  THE  LABEL  AND  CURRENT  VALUE  OF  A PARTICULAR  ITEM, ENTER 
A TRIPLE  (-RECORD  NUMBER, ITEM  NUMBER, 0) 

TO  PRINT  THE  ENTIRE  PARAMETERS  FILE  ENTER  THE  TRIPLE  -l.-l,0 
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WHEN  YOU  ARE  DONE  ENTERING  VALUES, ENTER  THE  TRIPLE  0,0,0 

>14,1,2239 

RECORD  14,  ITEM  1 

DAILY  CONTROL  MORTALITY  RATE^  FOR  PHASE  2 INSTAR  1 

HAS  THE  VALUE  .2239000 

>0,0,0 


INITIAL  CONDITIONS  FILE  GENERATION 

ENTER  NUMBER  OF  DAYS  STRESS  OPERATES  IN 

FIRST  INSTAR,  SECOND  INSTRAR,  THIRD  INSTAR,  FOURTH  INSTAR 

>0,0, 0,0 


TREE  CLASS  DEFINITIONS: 

YOU  HAVE  SPECIFIED  000  DOUGLAS  FIR  TREE  CLASSES 
AND  005  GRAND  FIR  TREE  CLASSES.  FOR  EACH  CLASS 
YOU  WILL  NEED  TO  ENTER  1 LINE  CONTAINING  THE  FOLLOWING  VALUES: 

1.  NOMINAL  % NEW  FOLIAGE 

2.  NOMINAL  TOTAL  FOLIAGE  BIOMASS 

3.  # OF  VIABLE  EGGS  PER  BRANCH 

4.  INSTAR  OF  EGG  COUNT 

5.  DAY  IN  INSTAR  OF  EGG  COUNT 

6.  ACTUAL  NEW  FOLIAGE  BIOMASS 

7.  ACTUAL  OLD  FOLIAGE  BIOMASS 


>30,230,50,1,1,69,161 

>30,230,75,1,1,69,161 

>30,230,110,1,1,69,161 


>30,230,150,1,1,69,161 
>30,  230,200,1,1,69,161 

IF  YOU  WANT  A LISTING  OF  THE  ENTIRE  *IC  FILE,  ENTER  -1,-1,0 

IF  YOU  WANT  TO  SEE  ONE  LINE,  ENTER  LINE  #,0,0 
IF  YOU  WANT  TO  CHANGE  A VALUE,  ENTER  LINE  #, ITEM#, NEW  VALUE 
TO  CONTINUE  WITH  PROGRAM,  ENTER  0,0,0 


>-l,-l,0 


2 

.00000 

.00000 

2 

. 00000 

. 00000 

2 

.00000 

.00000 

2 

. 00000 

.00000 

2 

.00000 

. 00000 

.0 

.0 

*IC.  AND  *PARAMETER. 
****GOODBYE**** 

>0FIN  9 


30.00000230.00000 

30.00000230.00000 

30.00000230.00000 

30.00000230.00000 

30.00000230.00000 


69.00000161.00000  50.00000 

69.00000161.00000  75.00000 

69.00000161.00000110.00000 

69.00000161 . 000001 50 . 00000 

69.00000161.00000200.00000 


FILE  PREPARATION  FOR  PHASE  2 START  COMPLETED 


1 

2 

3 

4 

5 


8 Direct  control  applied  at  Phase  II,  Instar  1. 

9 Cost  for  this  run  was  $1.30. 


17 


Finally,  run  a simulation  with  the  control  option  as  follows: 


Runstream  No.  5 


>0RUN  M01MAG, 1105203405  , CIESLA 

DUP  ID,  NEW  ID  IS  M01MAH 

DATE:  180680  TIME:  121251 

S2K  INFO  - , 52000,80-00006  (OINFO)MONDAY  11:35 

>0QUAL  WC2 

READY 

>0ADD,L  DFTM*RUNSTREAM.PRE2 


0DELETE,C  *LUN25. 

FURPUR  27R3A  E35  SL73R1  08/06/80  12:14:04 

LUN25  IS  NOT  CATALOGUED  OR  ASSIGNED 
FAC  STATUS:  400010000000 


0ASG,UP  *LUN  25. 
READY 


OUSE  25,*LUN25 
READY 


@DELETE,C  *LUN10. 

FURPUR  27R3A  E35  SL73R1  08/06/80  12:14:17 

@LUN10  IS  NOT  CATALOGUED  OR  ASSIGNED 
FAC  STATUS:  400010000000 


0ASG,UP  *LUN10. 
READY 


0USE  1O,LUN10 
READY 


0ASG,T  *4 
READY 


0ASG,T 

READY 


*5 


18 


PASG,T  *6 
READY 


0ASG,T  *7 
READY 


0ASG,T  *8 
READY 


0ASG,T  *9 
READY 


0ASG,T  *11 
READY 


PASG,T  *12 
READY 


0DELETE,C  *15. 

FURPUR  27R3A  E35  SL73R1  12:14:48 

*15  IS  NOT  CATALOGUED  OR  ASSIGNED 

FAC  STATUS:  400010000000 

@ASG,UP  *15 
READY 


OFREE  *15. 
READY 


0ASG,A  *15. 
READY 


0DELETE,C  *16. 

FURPUR  27R3A  E35  SL73R1  08/06/80  12:15:15 

*16  IS  NOT  CATALOGUED  OR  ASSIGNED 

FAC  STATUS:  400010000000 


0ASG,UP  *16 
READY 
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0FREE  *16. 
READY 


0ASG,A  *16. 
READY 


0DELETE,C  *17. 

FURPUR  27R3A  E35  SL73R1  O?/06/80  12:15:37 

* 17  IS  NOT  CATALOGUED  OR  ASSIGNED 
FAC  STATUS:  400010000000 


0ASG,UP  *17 
READY 


0FREE  *17. 
READY 


0ASG,A  *17. 
READY 


0ASG,T  *18 
READY 


0ASG,T  *19 
READY 


0ASG,T  *20 
READY 


0DELETE,C  *21. 

FURPUR  27R3A  E35  SL73R1  08/06/80  12:15:58 

*21  IS  NOT  CATALOGUED  OR  ASSIGNED 

FAC  STATUS:  400010000000 


0ASG,UP  *21 
READY 


0FREE  *21 
READY 


@ASG,A  *21. 
READY 


0ASG,T  *22 
READY 


0ASG,T  *23 
READY 


0ASG,T  *24 
READY 


PASG.A  DFTM*DFTMDF. 

FACILITY  WARMING  000200W000 


@USE  13,DFTM*DFTMDF 
READY 


(3ASG,A  DFTM*DFTMGF. 

FACILITY  WARNING  000200000000 

0USE  14,DFTM*DFTMGF 
READY 


0ASG,A  *IC. 
READY 


0ASG,A  *PARAMETERS 
READY 


@USE  30,*IC 
READY 


0USE  60,*PARAMETERS 
READY 

>0XQT  DFTri*OUTBREAK. 

>ADD,L  DFTM*RUNSTREAM.POST  2 
ADDX 

DATA  IGNORED  - IN  CONTROL  MODE 
>AD 

0ADD,L  DFTM*RUMSTREAM.P0ST2 


(3FREE  *2 

FACILITY  WARNING  1(210000000000 


(9FREE  *3 

FACILITY  WARNING  100000000000 


0FREE  *4 
READY 


0FREE  *5 
READY 


0FREE  *6 
READY 


0FREE  *7 
READY 


0FREE  *8 
READY 


0FREE  *9 
READY 


0FREE  *10 
READY 


0FREE  *11 
READY 


0FREE  *12 
READY 


0FREE  *13 
READY 


0FREE  *14 
READY 


0FREE  *18 
READY 
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0FREE  *19 
READY 


0FREE  *20 
READY 


0FREE  *21 
READY 


0FREE  *22 
READY 


0FREE  *23 
READY 


0FREE  *24 
READY 


0FREE  *25 
READY 

>0SYM,U  *15.,,FCR104 
>0SYM,U  *16.,,FCR104 
>0SYM.U  *17.,,FCR104 
>FIN  10 


10  Cost  for  this  run  was  $3.46. 


Model  ouputs  for  control  vs  no  control  under  these  conditions  are 
compared  in  tables  1-6  for  each  infested  area.  After  these  model 
simulations  have  been  made,  option  2 of  the  front-end  program  was  used  to 
create  the  specification  file  to  make  a post-simulation  run.  The  runstream 
for  creating  the  specification  file  follows: 

Runstream  No.  6 


>@RUN  M01MAG, 1105203405  ,W0NG 
DATE:  08/11/80  TIME:  115751 
>0ASG,A  DFTM*UTILITY1. 

READY 

>0XQT  DFTM*UTILITY1. 

DFTM  OUTBREAK  MODEL  DATA  PREPARATION  UTILITY 
VERSION  01.001  COMPILED  6/3/80 


PROGRAM  OPTIONS:  1=  IC  AND  PARAMETERS  FILE  GENERATION 
2=  SPECTAB4  FILE  GENERATION. 

3=  CONTROL  MORTALITY  ALGORITHM 
ENTER  OPTION  DESIRED 
>2 

*SPECTAB4.  FILE  GENERATION 
ENTER  RUN  QUALIFIER 
>WC2 

RUN  QUALIFIER  IS  NOW  WC2 
ENTER  NUMBER  OF  TREE  CLASSES 
>5 

ENTER  PAIRS  OF  NUMBERS(ONE  PAIR  PER  LINE)  FOR  EACH  ITEM  DESIRED 
FROM  THE  TABLE4.  PROGRAM. 

WHEN  FINISHED  ENTERING  PAIRS, ENTER  0,0 
>1,3 
>2.3 
>3.3 
>4.3 
>5.3 
>6.3 
>7.3 
>8.3 
>9.3 
>10.3 
>0,0 

OPTION  CODES: 

T=TERMIMATE  PROGRAM  NORMALLY 
C=CHANGE  AN  ITEM  IN  *SPECTAB4.  FILE 
A=ADD  AN  ITEM  TO  *SPECTAB4.  FILE 
L=LIST  CURRENT  *SPECTAB4.  FILE 
D=DELETE  AN  ITEM  FROM  *SPECTAB4.  FILE 
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ENTER  OPTION  CODE 
>T 

***WC2*SPECTAB4  HAS  BEEN  GENERATED  WITH 
*SPECTAB4.  FILE  PREPARATION  COMPLETED 
****GOODBYE**** 

>ofin11 


11  Cost  for  the  run  was  $0.87. 


ITEMS*** 


With  this  specification  file,  a computer  run  was  made  to  retrieve 
information  on  the  numerical  changes  of  the  population  resulting  from 
control.  The  runstream  for  the  post-simulation  run  follows: 


Runstream  No.  7 


>PRUN  M01MAG, 1105203405  ,W0NG 
DATE:  081180  TIME:  120252 
>0QUAL  WC2 
READY 

>(3ADD,L  DFTM*RUNSTREAM,TABLE4B 


0ASG,A  *SPECTAB4 
READY 


PASG.A  *LUN10. 
READY 


PASG,A  *LUN25. 
READY 


(3DELETE,C  *TABLE4. 

FURPUR  27R3A  E35  SL73R1  08/11/80  12:03:38 

*TABLE4  IS  NOT  CATALOGUED  OR  ASSIGNED 
FAC  STATUS:  400010000000 


0ASG,CP  *TABLE4. 
READY 


OUSE  11,*TABLE4 
READY 


OUSE  12,*SPECTAB4 
READY 


OUSE  10,*LUN10 
READY 


OUSE  25,*LUN25 
READY 
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0XQT  DFTM*NTABLE4.A2 

THERE  WERE  5 TREE  CLASSES  COUNTED  IN  *LUN10. 
REQUEST  IS  TO  PROCESS  THE  FIRST  5 TREE  CLASSES 
THERE  WERE  10  VALID  ITEMS  REQUESTED 


0FREE  11 
READY 


0FREE  *TABLE4. 

FACILITY  WARNING  100000000000 


PFREE  12 
READY 


0FREE  10 
READY 


0FREE  25 
READY 

>0SYM,U  *TABLE4.,2,FCR104 
>0FIN  12 


12  Cost  of  this  run  was  $0.53. 
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The  output  table  from  this  run  is  shown  in  Table  7. 

The  user  should  find  it  very  instructional  to  construct  examples  similar 
to  those  provided  in  this  report.  With  a few  practice  sessions  on  the 
terminal,  the  operating  characteristics  of  the  front-end  program  can  be 
grasped  quickly,  by  a person  with  very  little  background  in  computers.  When  a 
large  number  of  simulations  have  to  be  made,  this  program  should  prove  to  be  a 
real  time  saving  tool. 


PROGRAM  VARIABLE  DEFINITIONS 

In  this  Section,  the  definitions  of  the  variables  names  used  in  the 
program  source  codes  are  explained  in  terms  of  their  name,  type,  size,  and 
usage.  The  source  codes  for  this  program  can  be  made  available  upon  request. 
It  is  anticipated  that  additional  capabilities  will  be  installed  as  the  need 
arises.  Revisions  of  this  document  will  be  made  accordingly  to  provide  the 
user  with  the  latest  information. 


ICFILE  (Subroutine): 


NAME 

TYPE 

SIZE 

USAGE 

ITIC 

I 

200,9 

matrix  containing  the  integer  values  of  the  IC  file 

TIC 

R 

200,9 

matrix  containing  the  real  values  of  the  IC  file 

NOTE:  ITIC  and  TIC  have  been  equivalenced  and  are 
therefore  actually  stored  in  the  same  memory 
locations.  Use  of  ITIC  and  TIC  allow  both 
real  and  integer  values  to  be  stored  in  the 
same  matrix 

NDF 

I 

user-input  number  of  Douglas-fir  tree  classes 

NGF 

I 

user-input  number  of  Grand  fir  tree  classes 

NTC 

I 

total  number  of  tree  classes  (NDF  + NGF) 

NPHASE 

I 

user-input  starting  phase 

IND 

I 

counter  for  number  of  tree  classes  processed  so 
far.  In  this  version  also  use  as  the  tree  class 
index  (9th  item  in  the  current  record  of  the  IC 
file)  & equal  to  K 

TIC(K,8) 

R 

calculated  insect  density  for  Kth  tree  class.  If 
TIC(K,8)  is  greater  than  1000,  message  printed  and 
program  terminated 
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NAME  TYPE  SIZE 


USAGE 


ITIC(K,1)  I 

II  I 


12  I 


VAL  R 

IVAL  I 


species  code  (1  = Douglas-fir;  2 = Grand  fir) 

user- input  editing  code  for  IC  file.  May  be  a line 
(or  row,  or  record)  number  of  an  entry  in  the  IC 
file  or  merely  a code  for  printing  or  continuation 
of  program 

user-input  editing  code  for  IC  file.  May  be  an 
item  (or  column)  number  of  an  entry  in  the  IC  file 
or  merely  a code  for  printing  or  continuation  of 
program 

user-input  new  value  for  the  (Il,I2)th  entry  in  the 
IC  file  or  a code  for  printing  or  continuation  of 
program 

integer  form  of  VALUE  for  integer  entries  in  the  IC 
fil  e 


SDAY  I 6 Number  of  days  stress  operates  in  instars  1 to  6. 

SDAY  (1)  to  SDAY  (4)  are  user  input;  SDAY  (5)  and 
SDAY  (6)  are  always  zero. 

K I Line  (or  record)  number  of  current  tree  class  in  IC 

file 

PRFILE  (Subroutine): 


NTC 

NPHASE 

PRN 

IPRl 


I 


number  of  tree  classes 


I starting  phase 

R 24.6  matrix  containing  records  2-25  of  the  PARAMETERS 
file 


I 22  vector  containing  first  record  of  the  PARAMETERS 

fil  e 


ANSW 

ICODE 

I 

NPOSUB 


C*3 


character  variable  containing  a user-input  response 
(yes/no) 

numerical  code  corresponding  to  value  of  ANSW: 

"no"  = -1,  "yes"  = +1,  otherwise  = 0 

user-input  number  of  tree  class  subsets  before 
checking  for  allowable  values 

maximum  possible  number  of  subsets  given  the  number 
of  tree  classes  (NTC) 
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NAME  TYPE  SIZE  USAGE 


NSUBS  I 

number  of  tree  class  subsets  after  checking  for 
al  lowable  val  ues 

11  I 

current  user- input  low  subset  index  before  checking 
for  error.  Becomes  IPR1(J2  +1)  if  no  error 

12  I 

current  user-input  high  subset  index  before 
checking  for  error.  Becomes  IPR1(J2  + 2)  if  no 
error 

VALUE  R 

user-input  new  value  for  (11,  I2)th  entry  in  the 
PARAMETERS  file 

J I 

subset  currently  being  input  is  the  Jth  subset 

J2  I 

position  of  last  accepted  entry  in  IPRl  vector 
(=J*2  because  entries  1 and  2 in  IPRl  are  always 

MAIN  (Main  Program) : 

NPHASE  & NTC,  and  because  there  are  2 entries  per 
subset  pair) 

IPRl  I 22 

vector  containing  the  first  record  of  the 

PARAMETERS  file 

PRN  R 24,6 

matrix  containing  records  2-25  of  the  PARAMETERS 
fil  e 

lOPT  I 

option  code  for  desired  subroutine 

AMPLE  (Subroutine) 


I I 

record  (or  row)  number  of  item  in  PARAMETERS  file 

J I 

item  (or  column)  number  of  item  in  PARAMETERS  file 

IPRl  I 22 

vector  containing  first  record  of  the  PARAMETERS 
fil  e 

PRN  I 24,6 

matrix  containing  records  2-25  of  the  PARAMETERS 
file 

M I 

M = ABS(I) 

N I 

N = ABS(J) 

IR  I 

IR  = N module  2 (used  in  the  program  to  determine 
whether  an  item  number  is  odd  or  even) 
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NAME 

INDEX 

SALUTE 

IGEN 

IVER 

CDATE 

SOLONG 

lOPT 

NPHASE 

NCLASS 

NCL 

IL 

IH 

ABRTR 
IS  TAT 

GETFLS 

CCD 

QUAL 

QUALIF 


TYPE  SIZE 


USAGE 


I number  of  tree  class  subset  cooresponding  to  item 

(l,j)  in  the  PARAMETERS  file 

(Subroutine) : 

I generation  number  (changed  when  recompiling 

program) 

I version  number  (changed  when  recompiling  program) 

I 8 date  of  compilation  (changed  when  recompiling) 

(Subroutine) : 

I subroutine  option  code 

I starting  phase 

(Subroutine) : 

I number  of  tree  classes  (input  by  user) 

I lower  bound  to  number  of  tree  classes  (fixed  in 

program  version) 

I upper  bound  to  number  of  tree  classes  (fixed  in 

program  version) 


(Subroutine) : 

I error  status  code  returned  after  executive  request 

initiated  by  routine  FACSF 


(Subroutine) : 

C*40  contains  character  string  for  the  0QUAL  executive 

request 

C*72  user-entered  run  qualifier  including  leading  and 

trail i ng  bl anks 

C*1  12  qualifier  squeezed  down  to  12  characters  (including 

trailing  blanks,  if  any) 
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NAME 

TYPE 

SIZE 

USAGE 

QTEST 

C*1 

contains  sequential  characters  from  QUAE  for 
testing  (for  blanks).  Also  used  as  user-input 
response  code. 

IDUPMS 

0 

octal  error  status  code  for;  "another  file  already 
assigned  with  same  number  but  different  qualifer" 

I NOCAT 

0 

octal  error  status  code  for:  "no  cataloged  file 
with  name  given  in  an  ODELETE  executive  request" 

IS  TAT 

I 

error  status  returned  after  executive  request 
initiated  by  FACSF  routine 

EQUAL 

I 

counter  for  number  of  characters  in  qualifier 

ISTATR 

I 

same  usage  as  ISTAT 

QTEST 

I 

also  used  to  contain  a character- type  user  response 

I 

I 

counter  for  number  of  characters  in  QUAE  checked 
for  blanks 

TCLASS 

(Subroutine) : 

NEND 

I 

counter  for  number  of  EOF's  read  - if  greater  than 
2,  run  is  terminated 

NF 

I 

number  of  tree  classes 

TYPE 

R 

variable  indicates  record  number  of  appropriate 
stress  mortality  (DF  or  GF)  in  PRN  matrix--  = 18 
for  DF,  = 19  for  GF 

IPRl 

I 

22 

vector  containing  first  record  of  PARAMETERS  file 

PRN 

R 

24,6 

matrix  containing  records  2 through  25  of 

PARAMETERS  file 

NDFC 

I 

number  of  Douglas-fir  tree  classes 

FIC 

R 

200,9 

matrix  containing  entries  for  the  IC  file 

INSTAR 

I 

instar  of  insect  count 

DAY 

I 

day  within  instar  of  insect  count 

FIC6 

R 

if  NPHASE  is  greater  than  1,  6th  item  in  current 
record  of  FIC  (actual  new  foliage  biomass)  is  set  = 
FIC6 
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NAME 

FIC7 

I 

NPHASE 

D 

lENT 

ILEFT 

NDAY 

INITID 

CHANGE 

I 

J 

VAL 

IVAL 

IPRl 

PRN 

SPCT4 

IXl 

1X2 

I 


TYPE  SIZE 


R 


I 

I 

R 

I 

I 

I 6 

L 

(Subroutine) : 

I 

I 

R 

I 

I 22 

R 24,6 

(Subroutine) : 


I 


USAGE 

if  NPHASE  is  greater  than  1,  7th  item  in  current 
record  of  FIC  (actual  old  foliage  biomass)  is  set  = 
FIC7 

counter  for  number  of  times  user  has  input  a line 
of  data 

starting  phase 

total  survivorship  from  beginning  of  phase  to 
INSTAR,  DAY 

number  of  records  read  so  far 
number  of  records  left  to  be  read 
same  usage  as  SDAY  in  ICFILE  routine 
not  used  in  this  version 

indicates  record  (or  row)  number  of  entry  in 
PARAMETERS  file 

indicates  item  (or  column)  number  of  entry  in 
PARAMETERS  file 

new  value  for  (I,J)th  entry  in  PARAMETERS  file 
new  value  for  (l,j)th  entry  in  PARAMETERS  file 
vector  containing  first  record  of  PARAMETERS  file 
matrix  containing  records  2-25  of  PARAMETERS  file 

user-input  first  member  of  pair  before  checking  for 
error 


I 


user- input  second  member  of  pair  before  checking 
for  error 


I counter  for  number  of  characters  checked  for  blank 

in  QUAL.  Also  counter  for  number  of  pairs  entered 
in  II  & 12  vectors 

I number  of  tree  classes 


KTOTAL 
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NITEM 

I 

number  of  accepted  items  (pairs)  in  SPECTAB4 

NUM 

I 

position  number  of  pair  for  editing 

ICK 

I 

loop  index  for  checking  new  pair  for  duplication  in 
existing  pair  list 

CCD 

C*40 

contains  the  character  string  that  makes  up 
executive  request  "OQUAL  (qualifier)" 

QUAL 

C*72 

user-entered  run  qualifier  including  leading  and 
trailing  blanks 

QTEST 

C*1 

contains  successive  single  characters  from  QUAL  for 
testing  (for  blank  character).  Also  used  as 
user-input  response  code. 

QUALIF 

C*1 

12 

qualifier  squeezed  down  to  12  characters  (including 
trailing  blanks,  if  any) 

11 

I 

47 

vector  of  1st  members  of  SPECTAB4  pairs 

12 

I 

47 

vector  of  2nd  members  of  SPECTAB4  pairs 

IDUPMS 

0 

octal  error  status  code  for:  another  file  exists 
with  same  name  but  a different  qualifier  (usually 
following  an  PASG  executive  request) 

I NOCAT 

0 

octal  error  status  code  for:  no  existing  cataloged 
file  v;ith  name  given  in  an  0DELELTE  executive 
request 

IS  TAT 

I 

error  status  returned  after  executive  request 
command  initiated  through  FACSF  routine 

EQUAL 

I 

counter  for  number  of  characters  in  qualifier 
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